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Abstract 

The conditions under which time-energy uncertainty relations de- 
rived by Deffner and Lutz [9] for time-dependent quantum systems 
minimize the time necessary to excite such systems from their ground 
state to excited states are examined. The generalized Margolus-Levitin 
and Mandelstam-Tamm inequalities are worked out for specific fermionic 
and bosonic systems. 
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1 Introduction 

Time is an essential concept for the realization of fast electronic devices. As 
an example it may be useful to control and minimize the time needed by a 
two-level quantum system to jump from its ground state to an excited state. 
This time is constrained by Heisenberg's time-energy uncertainty relation. 
The optimization of the time duration of quantum jumps has been the object 
of numerous theoretical studies starting with Mandelstam and Tamm pQ and 
pursued up to present time [2j [3j HJ [5j EJ [7J |8] . 

The derived expressions concern stationary systems. Recently Deffher 
and Lutz [9] extended the Mandelstam- Tamm(MT) and the Margolus-Levitin 
(ML) inequalities to time dependent systems which are either intrinsically 
time-dependent or driven by an external time-dependent perturbation. 

In practice it is of interest to test the time needed by a quantum system 
to be driven from its ground state into excited states. In the present work 
the inequalities are analyzed and applied to physical systems to which their 
application may be of some practical utility. In the first part of the present 
work we examine the general conditions under which the ML inequality is 
optimized i.e. comes closest to an equality for specific values of the time 
interval over which the system evolves. In a second step the ML and MT 
relations are applied to fermionic and bosonic systems. 

In section 2 the inequalities are explicitly formulated and the conditions 
under which the inequalities can be optimized in the ML case are presented 
in section 3. Section 4.1 is devoted to a quantitative study of the TM and LM 
expressions which are applied to a fermionic Id quantum chain. In section 
4.2 the ML inequality is applied to a simple bosonic system coupled to an 
external time-dependent perturbation which can act in a weak or strong 
coupling regime. 
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2 Time-energy inequalities 

The general structure of the inequality can be written as |9J 



where Qo T = |(\E'(0)|\E'(r))| is the overlap between the wave function at time 
t — and t = r, Cor = arccosfior and AE(r) characterizes the energy- 
difference acquired by the system between the initial and final time. In the 
following h — 1. 

In the Mandelstam-Tamm formulation the energy denominator AE(t) of 
Equation (1) is given in terms of the variance of the energy 



AE(t) = 1 t dt l - /T / m t / xv — 2 

with 



<F(^>=^(t)|F(t)^(t)) (3) 

In the Margolus-Levitin formulation AE{r) =< E(t) > —E(0) where 
E(0) is the energy at t — and E(r) is the average energy of the system 
over a time interval [0, r] < -E(r) >= 1/r J E(t)dt. 



3 Extrema in the Margolus-Levitin expres- 
sion 

In the case of a time- independent system the realization of eq.(l) as an strict 
equality can be of great practical interest since it leads to the determination 
of the minimum time needed by the system starting from its ground state to 
excited states. This question has been successfully answered in [7]. 

The question may also be raised in the case where the hamiltonian dy- 
namics are time-dependent when it is of interest to find out the minimal 
time interval [0, r] needed in order to realize the equality. It comes out that 



3 



an analytic solution has not been found yet. A priori a less ambitious em- 
pirical answer might correspond to Cq t = and a minimum of the energy 
denominator AE(t). 

A rigorous answer to the optimization of the inequality towards an equal- 
ity consists in a determination of the maximum of R(t) which brings the 
r.h.s. of the expression closest or equal to r. 

The first derivative of this quantity with respect to r leads to 

w , h[AE(r)dC 0r /dr - C 0r dAE(r)/dr] 
dR{r)/dr = (4) 

An extremum is reached if dR(r)/dr = which leads to 
dCor/dr dAE(r)/dr 



C 0t AE(r) 
This extremum is a maximum if, for the corresponding value of r 



(5) 



AE(T)d 2 C 0T /dr 2 - Co T d 2 AE(r)/dT 2 < (6) 

A consequence of Eq.(5) can be observed if Cq t is maximized. This cor- 
responds to n 0r = |(*(0)|*(r))| = 0, the vectors |*(r)) and |^(0)) are 
orthogonal to each other. Then dC 0r )/dT = and induces dAE{r)/dr = 
if AE(t) ^ 0. But 



dAE(r)/dr=^l- <E ^ > =0 (7) 

T T 

which leads back to the expression of AE(r) if -E'(O) = 0. Hence the 
stationarity of AE{r) is correlated with the orthogonality of the vectors 
\ty(r)) and |\l/(0)). The stationary point r can be an inflexion point or an 
extremum. Then C$ T = n/2 mod(kn) and 

Whether or not this limit can be reached and > replaced by a strict equality 
depends on the system. 
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4 Models and applications 



4.1 Fermionic Id chain 

In a first step the time-energy inequality is applied to the time-dependent Id 
chain with even periodic boundary conditions already introduced in [TUl H] 

H = J/2(l + 7 )£ «. , + J/2(l - 7) £ ^f+i - to £ (7? (9) 
(0 (0 w 

where of is the x component of the Pauli matrix and similarly for the y 
and z components. The system is integrable and the wave function is given as 
a product of single particle wave functions with corresponding energies [T2] . 
The time dependence is generated by a local excitation of the last spin by 
an external magnetic field 

Hi N \t) = h 1 eM-t/r H )S z N (10) 

with Sff = a z N /2 pj. 

The wave function of the system is obtained perturbatively, up to second 
order in the interaction which works as a perturbation and leads to the 
expression of the overlap 

Q 0t = (0)|(1 + t/ (1) (0, r) + U^(0, r))|*(0)>| (11) 

where |^(0)) is the wave funtion at t = and U^(0,r), U^(0, t) are 
the first and second order contribution to the evolution operator U(0, r) = 
exp[—i(H T + C dtHi (t)]. A justification for the neglect of higher order 
contributions used in the numerical application is given in the Appendix. 

The average energy of the system over a time interval [0, r] reads 

E(r) = 1/r [ T m)\H(t)Mt)}/m)Mt)) - E(0) (12) 
^0 

where |\&(i)) is the perturbed wave function evaluated up to order 2 and E(0) 
the initial energy of the system. 
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(1) 


(2) 


R(r) 


0.1 


0.01 


0.16 


0.2 


0.01 


0.34 


0.5 


0.01 


0.04 



Table 1: R(t) as a function of 7 (l)and time th (2) 

4.1.1 General considerations concerning the application of the 
model 

The present external magnetic field on an unique spin state of the chain has 
been explicitly chosen to produce a weak effect on the chain in order to allow 
for a perturbative treatment. Under these conditions and for fixed J = 1 
it comes out that the wave function overlaps are not very sensitive to the 
strength of the magnetic fields , neither h nor h\ as long as these quantities 
stay in the range of unity. The same is true in the case of variations of 
the asymmetry parameter 7 which are fixed in the interval [0,1]. In all 
applications the length of the chain is N = 100. 

4.1.2 Application of the Tamm-Mandelstam expression 

Typical results concerning the r.h.s. of the TM expression are shown in 
Figure 1 as function of 7 with r = 100, h — 1, hi — 1. The results show 
that R{j) keeps very small when compared with r. 

4.1.3 Application of the Margolus-Levitin expression 

Using the ML expression leads to the results shown in Table 1 for R(r). Here 
r = 100, h — 1, hi — 1 in lines (1-2), h — 1, hi — 2 in line 3. R(r) keeps 
again quite small when compared with r. 

4.1.4 General comments 

As a consequence of the perturbative nature of the time dependence of the 
total Hamiltonian H(t) = H + H\ '(t) the overlap between the wave func- 
tions at t — and t — r staying in the numerator of R{r) keeps relatively 
sizable, hence Cq t never gets close to its maximum 7r/2 which would corre- 
spond to orthogonal states. 
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Figure 1: R(t) as a function of 7 for r=100. 
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The energy denominators in both the TM and ML expressions are at least 
an order of magnitude larger that the numerator so that R{r) gets a small 
quantity in every investigated case. Since the external excitation is confined 
to a small time interval the evolution of this quantity reaches quickly a quasi- 
stationary value after a very short number of time steps. The time-energy 
inequality is strictly verified. 

4.2 Bosonic model 

As a second application consider a simple time-dependent bosonic system 
with and 1 excitation quantum. The Hamiltonian is given by the expression 

H = A(a + a + 1/2) + V*(t)a + + V(t)a (13) 

where V(t) = Voe luJt is an external perturbation. The wave function which 
obeys the time-dependent Schroedinger equation is written as 

Mt)) = Y, c »® e ~ iEnt \ n> ( 14 ) 

n 

Here one keeps two states, the ground state \n >= |0 > and an oscillator 
quantum excitation a + \0 >. Then the normalized amplitudes c n (t) can be 
determined analytically as solutions of a system of two coupled first order 
differential equations 

dc (t)/dt = -ie- iu>t V(t) Cl {t)/h (15) 
dd(t)/dt = -ie luJt V*{t)c {t)/h 

The second order differential equation obtained for Co(i) shows three types 
of solutions depending on the sign of the determinant 

A = W 2 /h 2 - (w - A) 2 . 

The amplitudes c (t) and c\(i) show an oscillatory behaviour if A < 
or equal and a product of an exponential and oscillatory behaviour when 
A > 0. 

The knowledge of the time-dependent wave function 
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= c o (t)e- iEot \0 > +c 1 (t)e-^ 1 *a + |0 > 



(16) 



where Eq = A/2, E\ = 3 A/2 and co(t), c\(t) are the normalized ampli- 
tudes. The overlap Cq t and energy E{r) can be determined for any value of 



(*(0)|*(r)> = (i? (0,r) +i/ (0,r))e- <BoT + (^(0, r) + i/^O, r)) e - iJ5lT (17) 



where i?o(0, r) and Iq(0, t) are the real and imaginary parts of Cq(0)co(t) and 
similarly for i?i(0, r) and 7i(0,r). Then the expression of the energy reads 



£(r) = 1/r / A(|cx(t)| 2 + 1/2) + 2 J Re[\/*(t)c*(t)c (t)e l ^ 1 - £ »)Vt] - £(0)(18) 



4.2.1 General considerations concerning the model 

Here the strength of the coupling in the time-dependent part of the Hamil- 
tonian can be weak as well as strong: 

Three different regimes can be generated. They correspond to A > 0, 
A = 0, A < 0. 

4.2.2 Analysis of the Margolus-Levitin inequality in these differ- 
ent regimes 

Here the initial conditions are fixed as co(0)=l, cfco/dt=0. 

• If A < (weak coupling) one can come close to the maximum of C 0r 
which leads also closest to a state which is not far from being orthogo- 
nal to the initial state |^(0)}. Calculations indicate that R(t) reaches 
a maximum when Cq t gets to a maximum. This is in qualitative agree- 
ment with the considerations of section 3 about the fact that the de- 
nominator AE(t) of R(t) can reach an extremum when its numerator 
is maximum. 

Some values of R(r) are shown in Figure 2. Here A — 1, u — 2, 
V = 0.475. The oscillations are due to the oscillating nature of the 
external perturbation. 



r 




9 



1.17 


- 0(1.44) 


i 


weak coupling - 








6(122) 


R 








0.38 






0(0.47) - 


0.22 


i 


0(0.28) 
i 


i i 




50 


100 


150 200 



T 



Figure 2: R{r) as a function of r.The numerical values between parentheses 
in the figure correspond to C 0r (weak coupling case) 
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Figure 3: R(t) as a function of r where Cq t = 0.79 (strong coupling case) 
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• If A > (intermediate and strong coupling) the situation is similar to 
the one observed in the fermionic chain system. The final state \&(t) 
reaches a maximum of 7r/4, Cq t and R(r) reach very quickly a station- 
ary value. This can be explained by the structure of the wave function 
which stabilizes to a final value due to an exponential time component 
which multiplies an oscillatory term. 

Some values of R(t) are shown in Figure 3. Column 1 corresponds to 
the total time r, column 2 to the value of Cq t . Here A = 6, uj = 4, 
V = 3. 

In all cases the time-energy inequality is numerically verified. 



5 Conclusions 

The extension of the Margolus-Levitin and Mandelstam-Tamm time-energy 
inequality to time-dependent systems has been investigated on hand of two 
models, a fermionic Id chain and a two-state bosonic system. The numerical 
calculations are in agreement with the analytical expressions, a fixed time 
interval is always larger than the r.h.s. of the inequality expressions. 

In the case of the Margolus-Levitin formulation it has been shown that 
the minimum overlap between the wave function at the initial time and the 
time evolved wave function reaches a minimum when the r.h.s. of the inequal- 
ity is itself at an extremum. Numerical calculations confirm this point in the 
case where this extremum corresponds to an exact orthogonality or not. 



6 Appendix 

The perturbation induced by the external perturbation given in eq. (10) is 
essentially governed by the relaxation time t# and the field strength hi. To 
first order the strength of the perturbation is fixed by Tnh\. At second order 
the magnitude of the different contributions are governed by the strength 
factors 
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Si = h\r 2 H (19) 

5 2 = h\{l/D 2 ) 

5 3 = hl/(r H D) 
D~ (1/r^ + Ae 2 ) 

where Ae ~ eW — corresponding to single particle ground or ex- 

cited states. The single particle state energies are of the order of unity. In 
the numerical applications hi is chosen to be of the same order of magnitude 
and Tn is two to three orders of magnitude smaller. As a consequence the 
second order contributions are two to three orders of magnitude smaller than 
the zeroth order ones. This justifies the cut-off of the expansion at order 2. 
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